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The spectral functions are studied in conjunction with the dyadic Green’s functions for various 
media. The dyadic Green’s functions are found using the eigenfunction expansion method for 
homogeneous, inhomogeneous, periodic, lossless, lossy, and anisotropic media, guided by the Bloch- 
Floquet theorem. For the lossless media cases, the spectral functions can be directly related to the 
photon local density of states, and hence, to the electromagnetic energy density. For the lossy case, 
the spectral function can be related to the field correlation function. Because of these properties, one 
can derive properties for field correlations and the Langevin-source correlations without resorting 
to the fluctuation dissipation theorem. The results are corroborated by the fluctuation dissipation 
theorem. An expression for the local density of states for lossy, inhomogeneous, and dispersive 
media has also been suggested. 

PACS numbers: Valid PACS appear here 


I. INTRODUCTION 

The subject of thermal noise is of great interest. Ther¬ 
mal noise gives rise to Brownian motion, black body radi¬ 
ation, and noise in electrical circuit known as Langevin, 
Johnson, or Nyquist noise [1-4]. When a system is in 
thermal equilibrium with its environment, as much en¬ 
ergy is absorbed from the environment as is emitted to 
it. The fluctuation dissipation theorem (FDT) [5-8] both 
in classical and quantum form has been formulated to 
describe this physical phenomenon. Many recent works 
have expounded on this phenomenon. 

The recent interest in spontaneous emission, Casimir 
force (and the ability to measure it), has revived the in¬ 
terest in this area as well [9-21]. A closely related area 
is in near-held heat transfer where photons participate 
in energy transfer [22-26]. Because most electromagnetic 
environment is dissipative, and that quantum systems are 
assumed to be non-dissipative, the use of some formulas 
for Casimir force calculation has given rise to controver¬ 
sies: the most well known of which is the Lifshitz formula 
for Casimir force [11-13]. 

The subject of thermal equilibrium also emerges in 
quantum transport where electrons move from one reser¬ 
voir to another when the two reservoirs are not in equi¬ 
librium. This gives rise to the non-equilibrium Green’s 
function approach. Green’s functions and spectral func¬ 
tions are often dehned to describe the physics of quan¬ 
tum transport in micro- and nano-scopic devices [27, 28]. 
Even though spectral function is often used in quantum 
transport, its use in electromagnetics is less well known 
[29]. 


Lossless quantum systems are easier to study. The 
hrst quantization of electromagnetic helds were done in 
lossless systems, or weakly lossy systems [30-33]. The 
theory of quantum dissipation system has also been ac¬ 
tively researched [34-36]. These systems were studied 
with a quantum system in equilibrium with a thermal 
bath, sharing some commonality with FDT. Recently, 
FDT has been used to motivate the quantization of elec¬ 
tromagnetic field in lossy media [37-40]. 

In this paper, first, the Green’s dyadics [41, 42], and 
spectral functions are defined and derived in terms of 
modal expansions for electromagnetic system, ranging 
from simple homogeneous, lossless media to anisotropic, 
inhomogeneous, lossy media. The spectral function is re¬ 
lated to the local density of states (LDOS) [49, 50, 53] 
and the density of states (DOS). Physical interpretation 
is given to the electromagnetic spectral functions. Since 
the system is in time-harmonic oscillation, the energy dis¬ 
tribution of the harmonic oscillators can be related to the 
generalized Planck distribution. Finally, some equations 
are derived from LDOS and spectral function, which are 
also derivable from the FDT. 


II. HOMOGENEOUS UNBOUNDED MEDIUM 
DYADIC GREEN’S FUNCTION 

The dyadic Green’s function for a homogeneous 
medium is a solution to [41, 42] 

V X V X G(r, r') - k'^G{Y, r') = lJ(r - r') (1) 

where and w is the operating frequency. 

A corresponding vector wave function, which is also the 
eigenfunction, is defined to be the solution to 

V X V X F(k:, r) — k^¥{k, r) = 0 
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( 2 ) 



2 


with being the eigenvalue. The eigenvalues of this 
system are real because the operators are Hermitian or 
self-adjoint [42, 43]. For Cartesian coordinates, the vec¬ 
tor wave functions can be found easily. The general so¬ 
lution to (2) is of the form ae®^ *'. In an unbounded 
homogeneous medium, n can take on any value and is 
uncountably infinite in number. By the use of periodic 
boundary condition in a volume V = a x b x d, n can be 
discretized and made countably infinite in number. Then 

K = XKx + yHy + ZKz, 

Kx = ^iTTja, Ky = 2rmr/b, Kz = 2pT:ld, 

k'^ = K ■ K = ( 2 ? 7 r/a)^ -I- ( 2 m 7 r/ 6 )^ -I- ( 2 p 7 r/d)^ 

where {l,m,p) are integer triplets. To normalize the vec¬ 
tor wave functions, then 

H/T/ \ /rt\ 

M{K,r) = -(3) 

KsVV 

N(k, r) = —V X M(k, r) = ik X M(/«, r) (4) 

K 

= (5) 

where Kg = + Ky, and k = nf k. In the above, M is 

orthogonal to both z and k, N is orthogonal to M and 
K, while L is longitudinal and parallel to k. It can be 
shown that M, N, L above are orthonormalized, namely 
that 


I 5 dr — dll> dmm'^pp'^nn' (6) 

Jv 

where F„, n= 1,2,3 imply M, N, or L. 

In the above, k is an index associated with indices 
{l,m,p), and that M, N, L are mutually orthonormal. 
Hence, assuming the completeness of these eigenfunc¬ 
tions, then [55] 

Td(r—r') = 

[M(/«, r)a*(/«) -f N(/«, r)b*(/«) + L(k, r)c*(K)] 

K 

( 7 ) 


Using orthonormality, it follows that [56] 


horizontal, N as vertical (with respect to M and k), and 
L as longitudinal. By letting G(r, r') in (1) as 


K,n 


it follows that 


G(r,r',fc) 


E 


Fn{K,r)Fl{K,r') 

k-2 _ ^2 


( 10 ) 


( 11 ) 


where = k ■ k when n = 1,2, and K 3 = 0. Here, 
n = 1, 2 are used to index the transverse modes M and N, 
while n = 3 indexes the longitudinal mode L. The above 
is undefined unless a small infinitesimal loss is associated 
with k. 

The above is very similar to the solution of a modal 
cavity driven by a time-harmonic source with operating 
frequency lo where k = Lojc. The M and N modes are 
dynamic modes with resonant wave number k„, while 
the L modes are static modes with zero resonant wave 
number. But in this case, the modes are generated using 
periodic boundary condition rather than a cavity wall. 
In the above, (3) to (5) represent travelling waves while 
if cavity modes were used, they would represent standing 
waves. 


III. THE DENSITY OF STATES 

The density of states is the number of states (DOS) 
or modes per unit energy interval. For each eigenmode, 
there corresponds a frequency w = kc. Since photons are 
described here, the relation between photon energy and 
wavenumber is = fiion = Hnc, or k = £t^/{hc), h is 
the Planck constant, and c is the velocity of light. 

Since 


^21tt ^2rmr ^2pTT 

K = X -h y— -h z—— 

a b d 


( 12 ) 


each mode occupies a volume of {2'kY/[ abd) = Stt^/U in 
K-space. The number of modes or states in a spherical 
shell between k and k + An is 


AN = 2x Attk^Ak/{%tt^/V) = X^k^AkV (13) 

TT^ 


I(j(r — r') = [M(k, r)Ml(/«, r')-|- 

K 

N(/«, r)Nl(/«, r') -f- L(/«, r)Ll(/«, r')] (8) 


where one assumes that there are two polarizations for 
every mode. Notice that if AN is normalized by V, it is 
independent of V. 

A succinct way to write the density of states is 


or more concisely, 

I(I(r-r') = y^F„(/«,r)Fl(K,r') (9) 

K,n 

where F„, n = 1,2,3 denote M, N, or L, respectively. 
If z is assumed to be vertical, one can think of M as 


D{E) =Y,5{E - £n,n) (14) 

K,n 

where £K,n = knc, n = 1,2 represent the two polar¬ 
izations, horizontal and vertical, and E = Huj. Here, 
= \k\’^ is the eigenvalue of ( 2 ), which are non-zero 
for M and N modes. Since /« is associated with indices 
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{l,m,p), there are many degeneracies in the eigenvalue 
in (2): all k vectors from the origin to the sphere 
dehned by |/«| = k are degenerate. Also, in the above, 
one only retains states or modes with non-zero resonant 
wavenumbers. The L modes are not counted in the above 
since they have zero resonant wavenumber in (2). 

Equation (14) is like a spectroscopy function showing 
spikes whenever E coincides with one of the eigenstates 
with energy £n^n- The number of states in an energy 
interval E and E + AE is given by 

pE-\-^E 

/ D{E)dE = AN, {E < £K,n <E + AE) (15) 

JE 


It is seen that (14) is in fact the density of states per 
unit energy interval because integrating it over an energy 
interval gives the number of states in that interval, as 
shall be shown below. 

In the limit when V oo, one can convert the 
summation in (14) into an integral by the substitution 
that Ek = V/(2Tr)^ E ^^xAKyAKz = V/{2Tr)^ fdn. 
Namely, 


D{E) = 2 X 


V 


dhiS{E — Hkc) 


V 

47r3 


EiTK^dK6{E — hnc) 


Vk^ 

tt'^E 


(16) 


where it is noted again here that E = hkc = hu. In 
the above, = fine has no angular dependence, and 
hence, f dK = AiTK^dK. Consequently 


AN = D{E)AE = 


Vk^ 


Ak 


(17) 


which is the same as (13) except that one replaces k in 
(13) with k here since E = hoj is the counter. The above 
is self consistent with (14). 


IV. SPECTRAL FUNCTION AND LOCAL 
DENSITY OF STATES 


The photon (electromagnetics) local density of states 
(LDOS) is dehned as 

D{r,E) = Y,\Fn{K,r)\HiE-8K,n) (18) 

K,,n 

where n = 1,2 represent the two polarizations, and n 
indicates the plane wave direction of the eigenfunction 
since they are traveling plane waves. The above has the 
property that 

D{E)= I drD{r,E) (19) 

Jv 

since 

[ |F„(/^,r)|2dr = l (20) 

Jv 


The above implies that the LDOS, when integrated over 
r, yields DOS. Therefore, LDOS has the unit of states 
per unit energy per unit volume. It also associates differ¬ 
ent probabilities of hnding a photon in different locations 
denote by r. 

The LDOS can be related to the dyadic Green’s func¬ 
tion as shall be shown next. To this end, a spectral func¬ 
tion is dehned as 


A(r, r') = i 


G(r,r') 


G^(r',r) 


( 21 ) 


The second Green’s operator is the adjoint of the hrst 
one. The Green’s dyadic is an integral operator acting on 
a functional Hilbert space made of 3 vectors. The adjoint 
of this Green’s dyadic operator is dehned to satisfy the 
following relationship 


drf 2 (r)- / c?r'G(r, r') • fi(r') = 


dr' / dr 


it 


G (r',r).f 2 (r) ■ 


( 22 ) 


In this manner, the transpose of a scalar integral operator 
is such that g*(r, r') = g(r', r). This is similar to that for 
matrices, where matrix elements of transpose matrices 

are related by A* = [A] ... Hence, one swaps the 
L J ij 1'* 

argument of the second Green’s dyadic in (21) to make 
it the adjoint of the hrst Green’s dyadic operator. To 
elaborate, from the above, one obtains 

G“(r',r) = G'''(r,r'), or G“(r, r') = G'^(r', r) (23) 


In this manner, the spectral function so dehned above is 
a self-adjoint operator. 

It can be shown that the above represents a Hermi- 
tian system. For a reciprocal Green’s function, it further 
has the property that G (r,r') = G(r',r). Hence, the 
spectral function can be written as 


A(r, r') = i j^G(r, r') — G* (r, r')J = —2Qto [G(r, r')] 

(24) 

From (11), it can shown that 

F;(/«,r)F(,(K,r') 

^ I- j - 2^ - 2 _ ^*2- 

_ ^ F„(K:,r)Fl(K,r') 

^2 _ f,*2 


where = k- k, and the last equality is obtained by let¬ 
ting K —>• —K and one notices that F„(—/«, r) = F* (k, r). 
Moreover, it also follows from that the eigenvalues of (2) 
is real and hence, is real. Therefore 


A(r, r') = —2‘Am [G(r, r')] = i | G(r, r') — G*(r, r') 

1 


= -2^F„(,^,r)Ft(K,r')S 

K,n 


.2 _ ^2 
a- 


( 26 ) 
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One can rewrite 


A. Photon Energy Density 


1 _ 1 / 1 1 

2Kn\Kn — k + fc 

If /c —>■ /c + i ?7 to have a small imaginary part to mean 
loss, then 





irj 

+ rf 


V 

(Kn - kY + if 


(28) 


A photon with frequency w has quantization energy 
given hy E = hw. It is the result of energy quantization 
of the quantum harmonic oscillator. If the quantum har¬ 
monic oscillator is in thermal equilibrium with another 
system, it can be shown that the average photon number 
in the photon number state of the oscillator is [44, 45] 

” ^ ^ru^KkBT) _ I 

Then the average energy of the photon state is given by 


It can be shown that 

lim 2 ^ 2 = ^'^( 2 :) (29) 

»7-5-o -f 

Comparing, and using the above in (26), it means that 
A(r, r') = -29m [G(r, r')] = -2 X! F) 

K,n 

■- - [5{Kn - k) + S{Kn + k)] (30) 

2htn 

Furthermore, S{ax) = a~^S{x), Kn = £K,nlfic, and k = 
E/hc, imply that 

( ^2 ip ) = ^ 

(31) 


Consequently, 

A(r,r') = —29m [G(r,r')] 

= -TT E 

K,n 

Here, S{E + £i<i,n) is ignored since E -\- £i<i,n Y Oj and 
£K,n > 0 for all K. Another point to note regards the 
longitudinal modes (n = 3) with zero resonant frequency 
£K,n = 0. If one keeps track of modes only with £n,n > 0, 
they need not be included in the above expression. 

If r = r', and noticing that by taking the trace, 

one has Tr [F„(k, r)Ft (/t, r)] = X) F„i(/t, r)Fb(/t, r) = 

i—1 

|F„(K,r)p, it means that 

^2 _ 2^2 _ 

D{r, E) = -;^Tr [A(r, r)] = ;^Tr {9m [G(r, r)] } 

(33) 

where D{r,E) is the local density of states as given by 
(18). Hence, the LDOS can be found once the imaginary 
part of the dyadic Green’s function is known. In the 
above, the spectral function and the Green’s dyadic are 
implicit functions of to and hence, E. 


Q{E = hw)= (n+^^hw=^tuj+ _ i 

(35) 

The above result can be derived from the quantum sta¬ 
tistical mechanics of a quantum harmonic oscillator. In 
the limit when /i —>■ 0 , the above becomes the classical 
result showing that 0{E) = ksT: this is in agreement 
with the equipartition of energy theorem with O.^ksT 
per degree of freedom. A harmonic oscillator has two 
degrees of freedom, one for storing kinetic energy while 
another one for storing potential energy [46]. 

If T —>• 0, then Q{E) = (l/2)feu, which is also a result 
from quantum mechanics. It implies that the quantum 
harmonic oscillator has nonzero energy even if it is in the 
ground state: this is called the zero point energy. This 
zero point energy is there even when n = 0 , or with no 
photon [44, 45] is present on average. This is the source 
of vacuum fluctuation of the electromagnetic field, and is 
the source of the Casimir force at zero temperature [13]. 

Hence, if the density of states of the system is known, 
the average photon energy of the system is given by 

pOO 

{£) = / D{E)&{E)dE (36) 

Jo 

In the above, D{E)Q{E)dE is the average photon energy 
in the energy interval dE, and hence, when integrated 
yields the average energy of the system. If one lets E = 
tkj, and writes 


D{E)Q{E)dE = D{E)e{E)hduj (37) 


then one can identify, after using (16) for D{E) that 


U{cj) = —Q{E)h= — 


tt^E 


-huj 


huj 


gfia;/(fcBT) _ 2^ 


(38) 

is just the black body energy per unit volume per unit 
radian frequency according to the generalized Planck’s 
radiation law, except that the zero-point energy is absent 
in the original Planck’s law. 

If the local density of states is known, then the aver¬ 
age photon energy per unit volume or the photon energy 
density at the location r is given by 


poo 

{£{r))= D{r,E)&{E)dE 
Jo 


( 39 ) 
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Given that D{r, E)Q{E) implies energy per unit energy 
per unit volume, it is prudent to define a field correlation 
function [57] 

C(r,r',a;) = {£{r,uj)£\r' ,uj)) = —‘^^^A{r,r',uj)Q{huj) 

TT 

= <^rn [G(r,r',a;)] Q{hoj) 

TT ■' 

. 

where E = Huj. In this manner, the LDOS can be directly 
related to this field correlation function as 

^(r, E)e{E)dE = D{r,E)e{E)Muj = Tr[eoC(r, r, w)]da; 

(41) 

In the above, C has the unit of field square per radian fre¬ 
quency. The factor of 1/2 is not needed here because the 
above represents the energy density in both the electric 
field and the magnetic field. For a homogeneous medium, 
they are equal to each other. With the above, we can de¬ 
fine a density of states per unit radian frequency, D(r, uj), 
so that 

Z)(r, Lo)0{huj)duj = D{r, E)Q{E)dE = D{v, E)Q{E)hdw 

(42) 


or 


D{r, (jj) = D{r, E)h = TrSw [G(r, r, w)] (43) 

The relation expounded in (40) can also be confirmed in 
Section VIII, and are commensurate with FDT as shall 
be shown later. 


B. Physical Interpretation—Lossless Case 

The spectral function is formed by the difference of 
two Green’s function: it describes the response of a point 
source in equilibrium with its environment. If the Green’s 
dyadic describes a point source field that is leaving the 
source point, the adjoint of the Green’s dyadic describes 
a field that is returning back to the point source. If the 
Green’s function satisfies causality in the time domain, 
it is also called the retarded Green’s function, or the for¬ 
ward Green’s function. The adjoint Green’s function is 
non-causal, and it is also called the backward Green’s fuc- 
tion, the advanced Green’s function, or a time reversed 
Green’s function in the time domain. A forward Green’s 
function fore propagates the field from a source point to 
a field point, while a backward Green’s function back 
propagates the field from the field point to the source 
point. For instance, in the homogeneous medium case, 
one Green’s function yields an outgoing field from the 
source point, while the other one absorbs an incoming 
field to the source point. The reactive fields [42] that 
represent stored energy in these two Green’s functions 
cancel each other, and only the radiative fields remain in 
the difference. 


In the lossless medium, when the frequency is off res¬ 
onance, Kn ^ k, the forward Green’s function and the 
reversed Green’s function cancel each other exactly, but 
when Kn = k, they do not cancel each other giving rise 
to spikes in the spectral function. It yields the density 
of states of the system in (14) and the local density of 
states in (18). When the spectral function is multiplied 
by 0(i?), the average photon energy per state, it yields a 
quantity proportional to the energy density. Therefore, 
we have suggestively written the spectral function in (40) 
so that it is related to the correlation function of the field 
at two points. This relationship also follows from Sec¬ 
tion VIII and also the fluctuation dissipation theorem. 
Hence, the spectral function motivates a result that can 
be drawn from the fluctuation dissipation theorem. 


V. LINE BROADENING AND LOSS 


Notice that (14) resembles a set of spectroscopic lines 
that are infinitely sharp because there is no loss in the 
system corresponding to infinite Q. If loss is introduced, 
such spectroscopic lines will be broadened due to finite Q. 
To see this, the delta function in (28) will be broadened 
if the loss is not infinitesimally small. 

It is to be noted that (11) is still valid even when k 
corresponds to a lossy medium or is a complex number. 
The generalized spectral function then becomes 


A(r, r') = —2 Qto [G(r, r')] = i G(r, r') — G*(r, r') 
= -2'^Fn{K,r)Fl{K,r')j{Kn,k) 


where 

j{Kn, k) = 9m 


1 


„2 _ U2 

1 


2Kn \Kn-k Kn + k 


(44) 


(45) 


It can be shown that in the limit when V oo, the 
summation becomes an integral, giving [58] 

^j{Kn,k) = 2x J dKj{K,k) = ^k' (46) 


where k = k' + ik”, k' and k” are real numbers, but the 
integral is independent of k", the loss. The multiplica¬ 
tion by 2 is necessary because there are two polarizations 
per mode. In the above, k = where for a lossy 

medium, /tq and cq are complex numbers. But by the 
Kramers-Kronig relation [47], they have to be functions 
of frequency uj: their real and imaginary parts are related 
by Hilbert transforms [42]. 

A generalized local density of states can be defined as 
D{r, E)^Y. r)l2r(£:,^.„, E) (47) 

Kj,n 
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where 


T{SK.n,E) = —--f{Kn,k) 
ttE 

with K„ = £t^ n/{kc), E = huj, and c is the velocity of 
light in vacuum. In the above, as defined in (11) 
remains real when k is complex due to material loss. In 
this manner, 


lim r{£n,n,E)=6{E-£n,n) 
fe"—>-0 


Then a generalized density of states becomes 

D{E) = Y,^iS^^„,E) (48) 

K,n 


The above has a nice physical meaning since it can be 
shown that 


D{E)=Y,^i£K^n,E) 

K,n 


2k'^ 

ttE 


Kr,n 


V 2k'^ 
ttE 


dK'y{K, k) 


k'^V 

tt^E 


(49) 


Therefore, the spectral function can still be related to the 
LDOS in the lossy case. The DOS in this case is related 
to k', the real part of the k where the 7 (k, k) function 
peaks on the real axis of k. If we choose V = 1, then the 
DOS is in terms of the number of states per unit volume 
per unit energy. We will show in Section VIII that this 
relationship exists even for a general lossy, anisotropic, 
inhomogeneous medium. 


A. Physical Interpretation—Lossy Case 

The above formula is similar to (16) excepted that k 
is now replaced by k', the location on the complex k 
plane where 7 (k, k) peaks. The T function, which is the 
normalized 7 function, behaves like a broadened delta 
function, and in the limit of no loss case, reverts back to 
a delta function. In the lossless case, each k is associated 
with modes with k = /c as in (14) and (18), but in the 
lossy case, there is a cluster of modes with k values as¬ 
sociated with the peak of the 7 function when n = k' as 
indicated in (45) and (46). 

In a lossy medium, the first term of the spectral func¬ 
tion describes a decaying field, but the second term, a 
back propagating field, describes a growing field. This 
can be thought of as a lossy system in equilibrium with a 
thermal bath. The loss in the system is accompanied by 
Langevin sources [4] induced by the thermal excitation 
of the environment. These sources produce a field that 
grows, instead of decays with distance. 

In either cases, the spectral function describes a system 
in thermal equilibrium. In the lossless case, the system 
is in thermal equilibrium with sources at infinity, while 
in the lossy case, the system is in equilibrium with the 


Langevin sources due to the loss in the medium. In both 
cases, the system is in time harmonic motion because of 
thermal equilibrium: as much energy is supplied to the 
system as it is lost, so that the harmonic oscillators in 
the vacuum, atoms, and molecules are in simple time- 
harmonic oscillation. (Notice that the spectral function 
is Hermitian because it describes an energy conserving 
system.) 

But in the lossy case, when the system is excited to be 
in a simple time-harmonic motion, a cluster of modes are 
excited in unison to be in time-harmonic motion. When 
these harmonic oscillators are quantized to be associated 
with a quantum or packet of energy hw, the energy is 
distributed over the cluster of modes with different 
and weighted with different coefficients according to 7 . 
These modes can be given a probabilistic interpretation 
as in quantum mechanics. Because of the above physical 
interpretation, one can use (36) to (39) to obtain the pho¬ 
ton energy and photon energy density of a lossy system 
as well. 


VI. THE PERIODIC STRUCTURE CASE 

The periodic structures can be assumed to be perfect 
electric conductors (PEC) [59]. In this case, the electric 
vector wave functions or eigenfunctions are defined to be 
the solutions to eigenvalue problem 

V X V X Fe,„(K:, r) - K^re,„(K,r) = 0 (50) 

where k is the same as defined as above (3). By choosing 
a periodic PEC (or PMC) structure, the above eigen¬ 
value problem is similar to ( 2 ) except for the boundary 
condition to be satisfied on the surface of the PEC (or 
PMC). The above problem is easily proved to be Hermi¬ 
tian or self-adjoint, and hence, the eigenvalues are real. 
The above eigenvalue problem, in accordance with Bloch- 
Floquet theorem, will yield eigen-solutions where [53] 

F,,„(k, r) = ^ ^ (51) 

G X 

and G’s are the reciprocal lattice vectors, and A G 
{1,2,3}, with Gq, Gq, Gq forming an orthogonal triad 
unit vectors, and Gq is parallel to /« -|- G. Hence, 
A G {1,2} correspond to transverse modes, while A = 3 
corresponds to the longitudinal mode. 

The above is solved with the boundary condition that 
h X Fe,„ = 0 on the PEC surface of the scatterer where 
h is the unit surface normal. In this case, the eigenfunc¬ 
tions and eigenvalues have to be found numerically. It 
is to be noted that «:^(/«) is the eigenvalue of the equa¬ 
tion (50) with different choices of the Bloch-Floquet wave 
vector K. Hence, it is a function of k. 

It can be shown that Fe^„(K:, r) can be orthonormalized 
such that 

/ F}„(k;, r) • Fe,„/(/«', r)dr = (Ik,«/(!„„/ (52) 

Jv 
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where n indicates the band of the solution for the same 
K as shall be explained. 

In the homogeneous medium case, the transverse 
modes are degenerate, but in a periodic structure consid¬ 
ered here, they need not be [53]. Furthermore, the higher 
frequency modes are generated in the Brillouin zone due 
to the aliasing or wrapping of high k modes into the pri¬ 
mary zone. For unlike the homogeneous medium case, 
where n € {1, 2}, for a fixed k, there could be many 
for each k value where n G {1, 2, 3,...}. 

One can also make k countably infinite by having the 
periodic structure nested within a larger periodic struc¬ 
ture, so that K is discretized by a larger periodic bound¬ 
ary condition. One can also have only one unit cell within 
the larger periodic boundary condition. In this case, 
the problem becomes a single region scattering problem 
within a larger periodic boundary condition. 

A similar expression to (II) for the periodic structure 
case can be derived for the electric dyadic Green’s func¬ 
tion such that 


Ge(r,r') = 


Fe,„(/t,r)Ft^„(K,r') 
K? — k'^ 


(53) 


Due to reciprocity of the dyadic Green’s function, namely 
that Gg(r, r') = Ge(r',r), it implies that Fe,n(—r) = 

The magnetic vector wave functions can be derived by 
taking the curl of the electric vector wave functions. To 
make them orthonormal [42], one defines 

Fm.«('«,r) = —^VxFe,„(K,r) (54) 

Kn{K) 

where m above does not indicate an integer but ‘mag¬ 
netic’. It satisfies (50) but with a different boundary 
condition that h x V x F^.n = 0 since V x F^.n = KnFe,n 
[60]. The magnetic dyadic Green’s function is derived to 
be 


Gm(r, r') = 


2 -^ k -2 _ 1.2 

K.n 


(55) 


There is a need to include both types of Green’s func¬ 
tion, and the reason will be obvious when one derives 
the LDOS. We define the electric spectral function to be 


Ae(r,r')='i Ge(r, r') - Gg(r', r) 


—2'^m [Ge(r, r')] 

(56) 


where £K,n = kKnC, and Kn{i^) which is a function of k. 
Therefore, 

TT ^ 

Ae(r, 2,E) = -—Y Fe.n('«, l')Fln('«: r')6{E - £K,n) 

(58) 


A similar expression can be obtained for the magnetic 
spectral function. From the above, one derives the LDOS 
to be 


D{r,E) =- 2 ^Tr[Ae(r,r)-HA„(r,r)] 


2 ^ 

K,n 


|Fe_„(/«, r)|^ -h |Fm,„(K,r)|^ 6{E - 


(59) 

The above integrates over r to become DOS. One can 
also express the LDOS directly in terms of the Green’s 
dyadics, namely, 


Dir, E) = ^Tr {S5TO [Ge(r, r) -h G„^(r, r)] } 

D{y,uj) = hD{r,uj) 

= ‘^^°^° Tr {Sto [Ge(r,r) -H G^(r,r)] } (60) 


The LDOS involves both the electric and magnetic vec¬ 
tor wave functions. This is because LDOS gives the prob¬ 
ability of detecting a photon in a complex structure and 
this probability is proportional to the strength of the elec¬ 
tromagnetic field. In a periodic structure where standing 
wave can occur, the strength of the electric field may not 
have the same distribution as the strength of the mag¬ 
netic field. The above can be integrated over r to obtain 
the DOS. 

At this point, it is expedient to define field correlation 
functions as 


Ce(r, r',uj) = {£{r,u;)£\r',uj)) 

= -^^Ae(r,r',u;)0(/ia;) 

TT 

= Qfji [Ge(r,r',w)] 0(Sw) 


(61) 


Cm(r, r',w) = {'Hir,uj)'H\r',uj)) 
wee 
TT 

2,tO€Q 


= Am(r,r',w)0(ftw) 

TT 


Qto Gm(l’, I■^w) 0(fiw) 


(62) 


Then the LDOS can be related to these correlation func¬ 
tions as 


In a similar manner as before, one defines the spectral 
function, and show that 


1 


^2 _ U2 


ttE 


[SiE — £K,n) — <5(A -|- £K,n)] 


( 57 ) 


D(r,w)0(Sw)dw 


= Tr 


ieoCe(r, r, w) -K ^fioCmir, r, w) 


duj 


(63) 


In a homogeneous medium case, eoCe(r, r) = fi^Cmir, r), 
and the above reduces to (41). 



In general, k„(k) = K„(fc, is dependent on the 

direction of k, or anisotropic. This is unlike the homo¬ 
geneous medium case, Kn{i^) = k, where it is much sim¬ 
pler or isotropic. Therefore, when loss is included, the 
prove of Equation (46), valid for homogeneous medium, 
is more difficult for normalization now. But a general 
spectral function with line-broadening, local density of 
states, and density of states equations can be assumed, 
as shall be confirmed later. Moreover, (61) and (62) can 
be conhrmed in Section VIII and are commensurate with 
FDT. 


The electric dyadic Green’s function, despite the change 
in the orthonormality condition expressed by (66), can 
be dehned as before in (53), and the electric spectral 
function can be similarly defined. 

Similarly, a magnetic dyadic Green’s function can be 
defined as 

V X V^(r) • V X Gm(r,r') - ■ Gm{r,r') 

= K(r — r') (68) 


VII. INHOMOGENEOUS MEDIUM CASE 
A. Lossless, Reciprocal Case 


with its magnetic vector wave function to be the solution 
to 


In this case, the relative permittivity and permeability 
tensors er(r) = e(r)/eo and = A‘(i’)/Mo, respec¬ 

tively. They are Hermitian tensors [54]. The electric 
dyadic Green’s function satisfies the equation 

V X 7I“^(r) • V X Ge(r, r') - k'^er{r) ■ Ge(r, r') 

= T(5(r —r') (64) 

An electric vector wave function is defined to be the so¬ 
lution to 


V X ^(r) • V X Fe,„(K,r) - K^erir) ■ Fe,„(K,r) = 0 

(65) 


The medium can be arranged as a periodic structure, 
and the derivation of the eigenfunction above is guided 
by the Bloch-Floquet theorem. The eigenvalue «;«(/«) 
is then a function of k. To this end, one defines the 
eigenfunctions Fe_„(K, r) similar to the periodic structure 
case as in (51). But the Bloch-Floquet modes are er(r) 
orthonormal since we can prove from (65) that 


/ drFl^,{K,r) ■ er(r) •Fe_„(K',r) = Snn'SK, 
Jv 

Using the above, it can be shown that 
T5(r - r') = ^ Fe.„(/«, r)F[: „(k, r') • e,.(r') 

= ^er(r) • Fe,„(K,r)F];_„(/«,r') 


K' 


K,,n 


( 66 ) 


= ^ e? (r) • Fe,„(K, r) (r') • Fe,„(K,r') 


it 


(67) 


V X ^(r) • V X Fm,„(K, r) - ■ Fm,„(K, r) = 0 

(69) 


and its corresponding magnetic spectral function. 
The LDOS is related to the spectral functions as 


i)(r, w) = hD{r, E) = 


27r 


e=(r)-Ae(r,r)-e^=(r) 



(70) 


The modification is necessary because of the new or¬ 
thonormality condition expressed by (66). Upon substi¬ 
tuting (53) and (55) into the above, invoking the prop¬ 
erty of the vector wave function that follows (53), and 
going through the manipulation as in Section IV, equa¬ 
tions (26) to (32), finally, the LDOS can be shown to be 
modified as 


E) = -Y^ [rln('«, !■) • e,.(r) • Fe,n(K, r) 


K,n 


+rL,n(«> r) • r)] 6{E - 8k, n) (71) 


By using the orthonormality condition (66), the above 
integrates over r to become the DOS. The LDOS can 
also be written directly in terms of the dyadic Green’s 
functions as 


Dir,uj) = hD{r,E) = 


(r) • Sm [Ge(r, r)] • e?: (r) -|- (r) • [Gm(i’, r)] • (r)| 


(72) 


TT 
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With this modification, the photon energy density can 
be calculated accordingly. Notice that the LDOS be¬ 
comes large when er(r) or is large. Since the 

Kramers-Kronig relation [47] implies that a frequency de¬ 
pendent or dispersive medium is also necessary lossy, we 
assume that the medium is frequency independent. 

In the above, we can similarly define correlation func¬ 
tions for the electric and magnetic fields as in (61) and 
(62). Then we can relate them to the LDOS as 


D{r,uj)<d{Huj)doj = Tr 


-eoer"(r)-Ce(r,r,a;)-e=(r) 


+ 2 ^oMr (r) • Cm(r,r,w) • (r) 


dto 


(73) 


B. Lossy, Reciprocal Case 


The equation for the electric dyadic Green’s function 
in this case is the same as before. The compendium equa¬ 
tion for the vector wave function can be defined as be¬ 
fore in (65) and (69), but now, 7I,,(r) and er(r) are non- 
Hermitian, but symmetric tensors to represent a general 
lossy, reciprocal medium. In this case, the auxiliary or 
transpose equations have to be defined and k„(/«) is not 
a real number anymore [42, 48]. One can assume that 
the lossy medium forms periodic islands immersed in the 
lossless background. Then one can still use the Bloch- 
Floquet theorem as guidance in deriving the eigensolu- 
tion. Again, one can also use periodic boundary condi¬ 
tion to make k countably infinite. 

The Bloch-Floquet modes are (r) orthonormal since 
it can be shown that 

/ drF)]„^(K, r) • er(r) • Fe„'(K',r) = (74) 

Jv 

where Fg^n,a is the vector wave function of an auxiliary 
problem where material properties er(r) and /ir(r) are 
replaced by their conjugate transpose, implying an active 
medium. It can be shown that Fe^n,a = F* „ in this case. 
It can be shown that the electric dyadic Green’s function 
can be expanded as 


Ge(r,r') = 


Fe,n(K,r)F^_„(/t,r') 


K,n 


kI — 


(75) 


The above assumes the completeness of the eigenmodes 
even when the medium is lossy. The above is also 
analogous to the complex symmetric generalized eigen¬ 
value matrix systems in linear algebra where the right 
and left eigenvectors are orthogonal to each other via a 
matrix. [61] 

As before, the electric spectral function is defined to 
be 


Ae(r, r')=i Ge(r, r') - Gg(r', r) =-2Qm [Ge(r, r' 


( 76 ) 


The above form is more complicated compared to the 
lossless inhomogeneous medium case. Here, is com¬ 
plex in (75), and we see that the above gives rise to line 
broadening as before. Moreover, at one given frequency 
w, a cluster of modes is excited. However, as shall be 
shown later, the spectral function can be related to the 
correlation function of the field, and hence, LDOS. 

The permittivity and permeability tensors are of the 
form 


e = e' + ie", fj, = fx' + ifj," (77) 


where both the real and imaginary parts of the above ten¬ 
sors are real symmetric for reciprocal media and hence 
are Hermitian. For dispersive media, one can define effec¬ 
tive permittivity and permeability tensors to be [4, 51, 52] 


Cfi = 


d(jje' 

dto 


_ duju' 


(78) 


where the above tensors are Hermitian. Since the above is 
derived using perturbation argument, it is only valid for 
the low loss case. Similar to the lossless inhomogeneous 
case, the generalized LDOS is defined as 


D{r,(jj) = hD{r, E) 

= ^Tr (r) ’ [Ge(r, r)] • e| (r) 

-beo/lj (r) • [G™ (r, r)] • (r) | (79) 


The above LDOS is not provable for a general inhomo¬ 
geneous lossy medium, because the line broadening is 
rather complicated here. But we will justify it in the 
next section. 

The Kramers-Kronig relations [47] requires that a lossy 
medium is also dispersive, and the lossless medium be 
non-dispersive. Hence, the above reduces to (72) in the 
lossless medium case. The above analysis remains valid 
even if er(r) or pr(j') functions of w. As a conse¬ 
quence, the eigenvalue K„(/t,a;) is a function of both k 
and w, so are the eigenfunctions F„(/«, w, r) of the electric 
and magnetic types. Equations (74) and (75) are valid 
for dispersive media as long as w is fixed. 


VIII. SPECTRAL FUNCTION AND FIELD 
CORRELATION 

For the lossless media, we have shown that the spec¬ 
tral functions is related to the local density of states, and 
hence, they can be related to the field correlation func¬ 
tions. This is because the self correlation of the field 
is related to local energy density, and hence, the local 
density of states. 

For the lossy media case, it is not clear that the spec¬ 
tral function is related to the local density of states, but 
we can relate the spectral function to field correlation 
functions directly, as we shall show in this section. In 
this manner, the spectral function can be related to local 
energy density and hence, the local density of states. 
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A. Statistical Average of Random Fields 


Equivalently, the Green’s dyadic operator is then 


Before proceeding with this section, one needs to define 
some properties of the correlations of random fields and 
sources. Imagine that random Langevin sources are gen¬ 
erating random fields due to thermal excitation of these 
sources by the environment. These random fields are sta¬ 
tionary time processes, namely the time correlations of 
these fields depend only on time delay between the fields, 
or t — t' [8, 24]. Then the frequency correlations between 
these fields depend only on d{uj — uj'), namely, they are 
uncorrelated in frequency. When the Langevin sources 
are in thermal equilibrium with its environment, and as 
shall be shown, the correlation of the fields at different 
locations r and r' are related via the spectral function. 

If the Fourier transform of the random held E(r,t) is 
dehned as 


/ OO 

dtE(r,t)e“* (80) 

-OO 

then the E(r, w) has the dimension l/w times held. The 
general form of the correlation of the electric helds at two 
locations r and r' is 

(E(r, w)E'l'(r',a;')) = S{uj — w')(£(r, w)£^(r',a;')) 

= S{uj - u}')Ceir,r',u}) (81) 

Similarly, the correlation function for the magnetic held 
is of the form 


g, = {V^ k^Sr) ^ (85) 

We assume that the relevant boundary condition or ra¬ 
diation condition is imposed so that the inverse of the 
operator above is uniquely dehned. Using these nota¬ 
tions, an interesting expression for the spectral function 
operator can be derived. It can be shown that 



x>“ - x> - [el - Sr^ 


= 29m (X>) — 2A:^9m (£r) 


( 86 ) 


The above E would be zero for a lossless medium for 
which 7Ir(r) and er(r) are Hermitian. Hence E repre¬ 
sents the loss of the system. In the above, the adjoint 
operator denoted by superscript a, is as dehned in (22) 
where the 3-vector and Hilbert space nature of the space 
is accounted for. 

Multiplying the above by Qe and ^ from the left and 
right, respectively, the electric spectral function operator 
is given by: 


Ae=l{Qe-9‘l)=geTgl (87) 


(H(r,a;)H^(r', w')) = 5{uj — a;')('H(r, u;)'H^(r', w')) 

= d(w - w')C™(r,r',w) ( 82 ) 

In the above Ce and Cm are proportional to energy den¬ 
sity per unit radian frequency. 

B. Connection of the Spectral Function to Medium 

Loss 

The random helds in a medium at thermal equilibrium 
are generated by Langevin current sources [38]. To this 
end, one rewrites electric dyadic Green’s function equa¬ 
tion in operator notation as: 

(V - k^Sr) Qe=T (83) 

where 

^ ^ [V X 7I-i(r)Vx] , Sr ^ [er(r)] (84) 

These are the operators (or general Hilbert space rep¬ 
resentation of the operators) originally dehned in coor¬ 
dinate space. For lossy media, D and Sr are non Her¬ 
mitian, and the expressions right of the arrows are their 
coordinate space representations. 


The above can also be written more suggestively so that 

{p - k'^Sr) X = E (88) 

It implies that the spectral function operator is gen¬ 
erated by distributed sources on the right-hand side that 
are related to the Langevin sources. The right-hand side 
of the above can be interpreted as induced sources due to 
the back propagation of held via g^ into the lossy part 
of the medium represented by E. 

However, the spectral function operator can also be 
dehned more simply as Ae = i{Ge — 9e)- As shall be 
shown, the spectral function is proportional to the held 
correlation function, and it appears that the helds are 
correlated by the difference of a forward Green’s func¬ 
tion and a backward Green’s function, making it appear 
like a radiating source occurring concurrently with an 
absorbing source. 

It is to be noted that in (87), the operators in the 3- 
vector functional Hilbert space, and the f sign implies 
the adjoint operator dehned in this inhnite dimensional 
space. With this caution, the coordinate representation 
of the spectral function operator Ae in (87) becomes 
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Ae(r,r') = ^(r) J dr"V x Gr„(r,r") • ■ [29to/x^ ^(r")] • V" x Gg(r',r")] 

— J dr"Ge(r, r")- [2fc^9TOer(r")] •Gg(r',r") (89) 


In the above, one has made use of that 

[pr^ir) ■ V X Ge(r,r')]‘ = e~^{r') ■ V' x G™(r',r) 

(90) 

The above can be derived using reciprocity theorem [62]. 
It is the generalization of (1.4.14b) in [42] to inhomoge¬ 


neous, anisotropic media. Furthermore, it can be shown 
that 

{p~^) pI = 9m (pj.) (91) 

Therefore, (89) can be rewritten as 


Ae(r, r') = -29mGe(r, r') = -e^ ^(r) J dr” V x Gr„(r, r") • [29m/x^(r")] • [/x^ ^(r")]^' x Ge(r',r") 

— J c?r"Ge(r, r") • [2fc^9mer(r")] •Gg(r',r") (92) 


C. Connection of the Field Correlations to 
Langevin Sources 


Similarly, one defines the electric field due to magnetic 
Langevin sources as 


Alternatively, (92) above can be related to the corre¬ 
lations of electric fields. To this end, one defines electric 
field due to the electric Langevin sources as 


Ee(r, w) = iw/xo J cir'Ge(r, r', w) • J(r', w) (93) 

The average of the outer product of electric fields due to 
the electric type sources at two different frequencies is 
then given by the field correlation 


Em(r, w) = ^ 


dr"VxG^(r,r",w)-M(r",w) (95) 


(Ee(r,a;)E];(r',a;')) = J dr” J 

•(J(r",a;)Jt(r"',.;0)-Gl( 


dr"'Ge(r,r",a;) 
r',r'",w') (94) 


The average of outer product of electric fields due to the 
magnetic type sources is then 


J 


(E„(r,w)E)„(r', w')) = Cr ^(r) J dra J drt,V x Gm(r, r^, w) • (M(ra, w)M^(rb, w')) • [/x^ ^(rj,)] ^ • Vf, x Gg(r', r^, w) 


(96) 


In the above, one has made use of (90), the fact that 
G (r,r') = G(r',r) from reciprocity, and that p^. is a 
symmetric tensor. Consequently, assuming that the field 


due to electric Langevin sources is uncorrelated to the 
field due to magnetic Langevin sources, we have 
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(E(r,a;)E^(r',u;')) = (Ee(r, u;)E];(r', w')) + (E„(r, u;)E]„(r', w')) = J 
■ Ge(r',rh,w') + e7^(r) ' J J dnV x Gm{r,ra,uj) ■ (M(r 


Notice that the above is of the same structural form as 
(92). 

D. Connection of Spectral Function to the Field 
Correlations 

By comparing (92) and (97), and in order for them 
to be proportional to each other, it implies that the 
Langevin currents have correlations of the form 

(J(r,a;)J^'(r',a;')) 

= S{uj — uj')6{r — r ')—-Qm [er(r)] 0(fiw) 

(98) 


dva J drbGe{r, ra,uj) ■ (J(ro, u;)J'*'(rb,w')) 

a, w)M'l'(rf„ w')) • [p~'^irb)]'^ ■ V& x Gl{r',rb,uj) (97) 


and 


(M(r, u;)M''' (r', w')) 

= S{uj — uj')d{r — [/^^(r)] 0(fia;) 

TT 

(99) 


The above implies that Langevin sources are uncorrelated 
in frequency and space. It also implies that the field 
correlation function in (97) becomes 


(E(r,a;)E^'(r', w')) = d(w - w') je,, ^(r) J dr"V x Gm(i',r",w) • [/Xr(r)] G{fuv) ■ ^(r")] 

V" X Gg(r',r",u;) + f dr"Ge(r, r", w) • [er(i’)] 0(fiw) • Gg(r', r", w') ^ (100) 


r 


The above implies that the right-hand side is proper- are related as indicate before in equations (40), (61), and 
tional to the spectral function times Conse- (62). Namely, 

quently, the field correlation and the spectral function 


(E(r, u;)El(r', w')) = 6{uj — uj'){£{r,uj)S\r', oj')) = 6{uj — a;')Ce(r, r',oj) = 6{uj — oj') -9m [Ge(r, r')] Q{hoj) 

( 101 ) 


Or that Ce(r, r',u;) is connected to the electric spectral 
functions as 

Ce(r,r',a;) = ^^^9m [Ge(r, r')] Q{huj) (102) 

The assertions in equations (98) and (99) are justified 
for the following reasons: 

• The spectral function consists of two bilinear terms, 
one of which is linearly proportional to the permit¬ 
tivity loss and the other one linearly proportional 


to the permeability loss. We can also derive (98) 
and (99) by turning off the loss of each kind respec¬ 
tively and equating (92) and (97) with the appro¬ 
priate proportionality constant. 

• The Green’s dyadic is undefined for a lossless 
medium unless infinitesimal small loss is intro¬ 
duced. Hence, one thinks of (97), and hence, (100) 
being valid for the case of infinitesimally small loss 
media as well. Therefore, (102) is valid for the loss¬ 
less case as well. 
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• For the lossless case, (39), (61) and (62) are confir¬ 
mations that the above derivation leading to (102) 
is correct. Furthermore, the spectral function can 
be related to the local density of states, and hence 
energy density after multiplying by 0(Sw). This is 
further affirmation that Equation (102) is at least 
valid when r = r'. 

A correlation function for the magnetic field can simi¬ 
larly be defined as 

(H(r, u;)H'l'(r',w')) = 5(w - w')Cm(r, r', w) (103) 

where Cm(r, r', uj) can be connected to the magnetic spec¬ 
tral functions as 

Cm(r,r',uj) = [Gm(r, r')] 0(/iw) (104) 

The above results are also corroborated by the fluctua¬ 
tion dissipation theorem [8, 24]. 

E. Physical Interpretation 

Now, however, (92) and (100) have a nice physical in¬ 
terpretation. It implies that there are two types of con¬ 
tributions to the electric field correlation function or the 
electric spectral function: one that comes from magnetic 
type Langevin current sources and the second that comes 
from electric type Langevin current sources. Since the 
complex conjugation of a field corresponds to time rever¬ 
sal or back propagation, a field at point r' back prop¬ 
agates into the lossy medium to the source point at r" 
via the backward Green’s function Gg. The Langevin 
source at point r" can only correlate or be coherent with 
the same Langevin source. The same Langevin source 
radiates a field that propagates from the source point r" 
to the field point at r via the forward Green’s function 
Gg. Different Green’s functions are used depending on 
the source type. The first term corresponds to magnetic 
type sources while the second term corresponds to elec¬ 
tric type sources. 

As can be seen, the above shows that the fields at two 
locations r and r' are correlated if they come from the 
same Langevin source. Most important, the above also 
shows that the two fields at two locations r and r' are 
correlated directly by the spectral function as in the left- 
hand side of (92), multiplied by &{Iiuj) and a multiplica¬ 
tive constant, which is a much simpler representation of 
the correlation function. The spectral function is non¬ 
singular at r = r' and is seen to be set up by distributed 
Langevin sources. 

F. Energy Correlation Functions 

In the above, we have derived the field correlation func¬ 
tions, and show their relationships to the spectral func¬ 
tions. To be more precise, when the energy density is 


needed, the energy correlation function needs to be de¬ 
rived as follows. An energy correlation function for the 
electric field, in connection to (78), is then defined as 

(eJ(r).E(r,w)Et(r',w')-el(r')) 

= (I(a; — u;')We(i’, I■^w) (105) 

In the above, We(i',i'',w) should have the dimension of 
energy density per radian frequency, and it can be con¬ 
nected to the electric spectral functions as 

We(r,r',a;) = (r) • 9m [Ge(r,r')] ■ei{r')Q{fuv) 

(106) 

The above can be related to the LDOS when r = r'. 

Similarly, a correlation function for the magnetic held 
is dehned as 

= S{u} - w')Wm(r,r',w) 

(107) 

With the understanding that WTn(r, r', w) should have 
the dimension of energy density per radian frequency, it 
can be connected to the magnetic spectral functions as 

Wm(r,r',w) 

= — 71 I (r) • 9m [Gm(r, r')] • 71 I (r')0(fiw) 

TT 

(108) 

The above can be combined and are related to the LDOS 
dehned in (79). 

IX. CONCLUSION 

The Green’s dyadics in terms of eigenmode expansions 
have been derived for general media including lossy, inho¬ 
mogeneous, anisotropic media. Then the electromagnetic 
spectral functions are dehned. The spectral function as 
dehned has the physical meaning of describing a source in 
equilibrium with its environment. It consists of a causal 
or radiative part, subtracting a non-causal, or absorptive 
part, representing the radiation and absorption of helds 
by Langevin sources. 

For lossless media, it can be shown that the spectral 
function is related to the local density of states. For lossy 
media, due to line broadening, the spectral function de¬ 
scribes a source at w that excites a cluster of modes in 
time harmonic motion, even though these modes may 
not be degenerate. Furthermore, by relating the spec¬ 
tral function to the held correlation function, the spectral 
function, including that for lossy anisotropic inhomoge¬ 
neous media, can be related to the local density of states. 

Many relations derived using FDT can also be derived 
using spectral function. One of them is that the corre¬ 
lation of the helds at two points is related via the spec¬ 
tral function which is the imaginary part of the Green’s 
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function. The second one is that the correlation of two 
Langevin currents are related via the imaginary part of 
their permittivity or permeability depending on the cur¬ 
rent type. Hence, our results are commensurate with 
FDT: namely, the same conclusions from spectral func¬ 
tions can be arrived at from FDT. It also justifies the 
use of FDT for the calculation of Casimir force in a lossy 
medium. 
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